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Abstract 


A differential detection technique for MPSK, which uses a multiple symbol 
observation interval, is presented and its performance analyzed and 
simulated. The technique makes use of maximum-likelihood sequence 
estimation of the transmitted phases rather than symbol-by-symbol detection 
as in conventional differential detection. As such the performance of this 
multiple symbol detection scheme fills the gap between conventional (two- 
symbol observation) differentially coherent detection of MPSK and ideal 
coherent of MPSK with differential encoding. The amount of improvement 
gained over conventional differential detection depends on the number of 
phases, M, and the number of additional symbol intervals added to the 
observation. What is particularly interesting is that substantial performance 
improvement can be obtained for only one or two additional symbol 
intervals of observation. The analysis and simulation results presented are 
for uncoded and trellis coded MPSK. 
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1.0 Introduction 

It is well known that, in applications where simplicity and robustness of 
implementation take precedence over achieving the best system performance, 
differential detection is an attractive alternative to coherent detection. Aside 
from implementation considerations, it is also possible that the transmission 
environment may be sufficiently degraded, e.g., a multipath fading channel, 
that acquiring and tracking a coherent demodulation reference signal are 
difficult if not impossible. Here again, differential detection is a possible, and 
perhaps the only, solution. 

In the past, differential detection of multiple-phase-shift-keying (MPSK) 
has been accomplished by comparing the received phase in a given symbol 
interval with that in the previous symbol interval and making a multilevel 
decision on the difference between these two phases [1]. An implementation 
of such a receiver and the analysis of its error rate performance on an additive 
white Gaussian noise (AWGN) channel may also be found in [2: Chap. 5]. In 
arriving at the results in [1,2], the assumption was made that the received 
carrier reference phase is constant over at least two symbol intervals and thus 
has no effect on the decision process when the above-mentioned phase 
difference is taken. This assumption is crucial to the analysis but is also 
realistic in many practical applications. Also, since the information is carried 
in the difference between adjacent received phases, the input information 
must be differentially encoded before transmission over the channel. 

Although differential detection eliminates the need for carrier acquisition 
and tracking in the receiver, it suffers from a performance penalty (additional 
required SNR at a given bit error rate) when compared with ideal (perfect 
carrier phase reference) coherent detection. The amount of this performance 
penalty increases with the number of phases, M, and is significant for M > 4. 

For example, at a bit error probability = 10'^, differentially detected BPSK 
(often abbreviated as DPSK) requires about 0.75 dB more bit energy-to-noise 
ratio (E^/Nq) than coherently detected BPSK. For QPSK (M = 4), the 
difference in E^/Nq between differential detection and ideal coherent 
detection at = 10'^ is about 2.2 dB. Finally for 8PSK, the corresponding 
difference in E^/Nq performance between the two is greater than 2.5 dB. 


Thus, it is natural to ask: Is there a way of enhancing the conventional 
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(two symbol observation) differential detection technique so as to recover a 
portion of the performance lost relative to that of coherent detection, and yet 
still maintain a simple and robust implementation? Furthemore, if this is 
possible, what is the tradeoff between the amount of performance recovered 
and the additional complexity added to the conventional differential 
detection implementation? The answers to these questions stem from the 
idea of allowing the observation interval over which symbol decisions are 
made to be longer than two symbol intervals while at the same time making 
a. joint decision on several symbols simultaneously as opposed to symbol-by- 
symbol detection. As such, one must extend the previous assumption on the 
duration of time over which the carrier phase is constant to be commensurate 
with the extended observation interval. For observations on the order of 
three or four symbol intervals, this is still a reasonable assumption in many 
applications. 

The theoretical framework in which we shall develop this so-called 
multiple-bit differential detection technique is the maximum-likelihood 
approach to statistical detection. In the next section, we derive the 
appropriate maximum-likelihood algorithm for differential detection of 
uncoded MPSK and show how the conventional technique is a special case of 
this more general model. Since, as mentioned above, we will be making joint 
symbol decisions in this new configuration, the technique is a form of 
maximum-likelihood sequence estimation, although no coding of the input 
information is implied. Later on in the report, we extend the theory 
developed here to the case of trellis-coded MDPSK. 

2.0 Maximum-Likelihood Detection of MPSK over an AWGN Channel 

Consider the transmission of MPSK signals over an AWGN channel. 

The transmitted signal in the interval kT < t < (k+l)T has the complex form 

s k = yf2Pe J * k (1) 

where P denotes the constant signal power, T denotes the MPSK symbol 
interval, and <j>k the transmitted phase which takes on one of M uniformly 
distributed values (i m = 27tm/M; m = 0, 1, ..., M-l around the unit circle. The 
corresponding received signal is then 

r k = s k e j0k + n k (2) 

where n^ is a sample of zero mean complex Gaussian noise with variance 
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and 0^ is an arbitrary phase introduced by the channel which, in the absence 

of any side information, is assumed to be uniformly distributed in the 
interval (-7 C,tc). 

Consider now a received sequence of length N and assume that 0^ is 
independent of k over the length of this sequence, i.e., 0^ = 0. Analogous to 
(2), the received sequence r is expressed as 


r = §e l9 + n (4) 

where r^, sfc, and n^ are, respectively, the kth components of the N-length 
sequences r, s, and n. For the assumed AWGN model, the a posteriori 
probability of r given £ and 0 is 


P(r\s,6) = - — j 

(2 nal) [ 2 < t „ 


\\r-se‘ 


i=0 

Simplifying the right hand side of (6) results in 


(5) 

( 6 ) 


where 



(7) 

( 8 ) 
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Since 0 has been assumed to be uniformly distributed, then the a 
posteriori probabability of z given s is simply 


p(rls) = J x P(r\s,0)p(6)dd 


1 


7- 2vW? eX PI 




\N-l 




1=0 


(9) 


where Iq(x) is the modified Bessel function of the first kind. Note that for 
MPSK, I sfc I ^ is constant for all phases. Thus, since I(j(x) is a monotonic 
function of its argument, maximizing p(r I $) over £ is equivalent to finding 


w-i 


max 

X 


-iCi 
<= o 


which, using (1), results in the decision rule 


(10) 


choose 0 if 


N - 1 

I 

1=0 






is maximum 


(id 


* 

where ^ is a particular sequence of the P m ’s. Note that this decision rule has 

a phase ambiguity associated with it since the addition of an arbitrary fixed 

phase, say <(> a , to all N estimated phases <j> k ,<p k _ l ,-..,<j> i _ N + l results in the same 

decision for <fc. Thus, letting <|> a = <J>k-N+l' the above decision rule can be 

* 

alternately expressed as choosing the sequence </> that maximizes the statistic 


;=o 


( 12 ) 


To resolve the above phase ambiguity, one should differentially encode 
the phase information at the transmitter. Letting 


4> k = + A$ k (13) 

where now A<j>k denotes the input data phase corresponding to the k th 
transmission interval and <f>k the differentially encoded version of it, then 
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N-i-2 

Qk-i ~ Qk-N+l = ^jAQk-i-m 
»»= 0 

and the above decision statistic becomes 


(14) 


V = 


N-i-2 


N-2 




r k-N + 1 + 






m=0 


i=0 


(15) 


This statistic implies that we observe the received signal over N symbol time 
intervals and from this observation make a simultaneous decision on N-l 
data phases. 


Some special cases of (15) are of interest. For N = 1, i.e., an observation of 
the received signal over one symbol interval, (15) simplifies to 


V=\r k f < 16 ) 

which is completely independent of the input data phases and thus cannot be 
used for making decisions on differentially encoded MPSK modulation. In 
fact, the statistic of (16) corresponds to the classical case of noncoherent 
detection which is not applicable to phase modulation. 


Next, let N = 2, in which case (15) becomes 


n = H-, + = k-,f + + 2Re{ Vl >-*»‘} (17) 

This results in the well-known decision rule for conventional MDPSK, 
namely, 

choose A<j> k if Re|r Jk r k *_ 1 e”^**J is maximum (18) 

which is implemented in complex form as in Figure 1. Thus, we see from 
this approach that conventional differential detection of MPSK is the 
optimum receiver in the sense of minimizing the symbol error probability 
given that the unknown carrier phase is constant over two symbol times. 

This result is not new other than, perhaps, the approach taken to demonstrate 
it. 


Now, to see a new structure, we consider (15) for N = 3. Here we have 
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77 = \r k _ 2 + r k e- J(A * k+Mk - l) + r k _,e )A * k ~ f = 

U 2 + ki-il 2 + tf + 2 Re{r 4 r 4 V-^‘^‘-»} 

+ 2 Re{r 4 _, r k _ 2 e~ JA * k ~ 1 ^ + 2 R (19) 

Thus, the decision rule becomes 

choose A<j> k and A$ k _ x if Re { r * r *- i e JA * k + r k-i r k- 2 e ~ jA * k l * + r k r k- 2 e /( ^* + ^*- l) | is maximum 

( 20 ) 

Note that the first and second terms of the metric used in the decision rule of 
(20) are identical to those used to make successive and independent decisions 
on A <p k and A<f> k _ x , respectively, in conventional MDPSK. The third term in 
the optimum metric is a combination of the first two and is required to make 
an optimum joint decision on A <p k and A 0 t _,. 

Clearly, a receiver implemented on the basis of (20) will outperform 
conventional MDPSK. Before demonstrating the amount of this performance 
improvement as a function of the number of phases, M, we first discuss the 
implementation of the optimum N = 3 receiver. Figure 2 is a parallel 
implementation of the decision rule of (20). It should be noted that the m 2 
phasors 1 needed to perform the phase rotations of the output r^r^)* can 
be obtained using a matrix which performs all possible multiplications of the 
M phasors e"jPo, e'jPl, ..., e'jPM-1 with themselves. Figure 3 is a series 
implementation of the same decision rule which, although simpler in 
appearance than Figure 2, requires envelope normalization and additional 
delay elements. 

3.0 Bit Error Probability Performance 

To obtain a simple upper bound on the average bit error probability, P^, of 
the proposed N-bit detection scheme, we use a union bound analogous to that 
used for upper bounding the performance of error correction coded systems. 

In particular, the upper bound on P^ is the sum of the pairwise error 
probabilities associated with each (N-l)-bit error sequence. Each pairwise 
error probability is then either evaluated directly or itself upper bounded. 

1 1n reality, only M phasors are needed since the sum angle when taken 

modulo 2 jc ranges over the set {Jq, Pi,—, Pm-1* 
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Mathematically speaking, let A0 = (A0 4 ,A0 jk _ 1 ,...,A0 4 _ w+2 ) denote the sequence 

A A A A 

of N-l information phases and A# = (A0 t ,A0 t _, t ...,A0 4 _ w+2 ) be the 
corresponding sequence of detected phases. Let u be the sequence of b = 
(N-l)log2M information bits that produces A0 at the transmitter and u the 

sequence of b bits that result from the detection of A$. Then, 

(21) 

where w(u,u) denotes the Hamming distance between wand u and 
Prjfj > rj I denotes the pairwise probability that A^ is incorrectly chosen 

when indeed A</> was sent. The decision statistic T] is defined in (15) and the 
corresponding error statistic T] is identical to (15) with each A <j> k replaced by 
A<j) k . For symmetric signalling sets (such as MPSK), (21) satisfies a uniform 
error probability (UEP) criterion, i.e., the probability of error is independent of 
which input phase sequence A0 is chosen as the correct sequence. Under 
these conditions, (21) simplifies to 

P b * 717 ui 77 X w ^>^) Pr {n > C 22 ) 

(W-l)log 2 M 1 — J 


where A<p is any input sequence (e.g., the null sequence (0,0,...,0) = 0) . 

3.1 Evaluation of the Pairwise Error Probability 

To compute Prjfj > rjIA^J, we use the approach taken in [3] for evaluating 
the performance of noncoherent FSK. It is convenient to define 


N-i- 2 

N-i ^ N-l 

w(A<j>) = X r t-. g " = ° I w(A(j)) = Y, r k-, e " =0 


(23) 


i=0 


1=0 


in which case. 


77 = |iv(A^)| 2 ; f) = w(A0 )| 2 


(24) 
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Then, the pairwise error probability Prj fj > i)\ Jis derived in Appendix A as 

Pr{fj>Ti\Ap} = ^[l-Q(4b,^) + Q(^,^b)] (25) 

where Q(a,(3) is Marcum's Q-function [4] and 





with E s = PT denoting the energy per data symbol and 


(26) 


. W-l ‘ 

<5= 


N-i-2 


A W-l * ^ 
= 2> 


(27) 


i=0 i=0 

In (27), it is understood that the summation in the exponent evaluates to 
zero if the upper index is negative. 


Note that for any given N, M, and input data sequence A0, 5 can be 
evaluated for each error sequence A0 . We now consider the evaluation of 
(22) and (25) for some special cases. 


3.2 Case 1: Conventional DPSK (N = 2, M = 2) 


From (27), we immediately get 5 = 0 and thus from (26) 



Substituting (28) into (25) gives 


(28) 


Pr{77>77ld0} = ^ 


1 ~Q 


I2E. 


N, 


,0 


o J 


+JII 


From the definition of the Q-function, 


(29) 


Q(a, 0) = 1; 


<2(0,0) = exp 



(30) 


Since for the binary case the pairwise error probability is indeed equal to the 
bit error probability, we have from (29) and (30) that 
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1 ( F ^ 

P b = —exp — *- (31) 

* 2 N 0 ) 

which is the well-known result for DPSK. 

3.3 Case 2: N = 3. M = 2 

Here there are three possible error sequences of length 2. The pertinent 
results related to the evaluation of (26) and (27) are given below: 



Since the Hamming distance vv(u,«) is equal to 1 for the first two error 
sequences and is equal to 2 for the third sequence, then using (32) in (25) and 
(26), the upper bound on bit error probability as given by (22) is evaluated as 



(33) 


To see how much performance is gained by extending the observation 
interval for differential detection from N = 2 (conventional) to N = 3, we 
must compare (33) to (31). Due to the complex form of (33) this comparison is 
not readily obvious without resorting to numerical evaluation. On the other 
hand, by examining the asymptotic (large E s /Nq) behavior of (33) we can get 
an immediate fix on this gain. 

When both arguments of the Q-function are large, the following 
asymptotic approximations are valid [6]: 
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P_ 
P V2^a 


exp 


ioc-1 y 


Q(a,p) = 


1 


/? - a V 2na 


exps 


2 

V 

( P ~ «) 2 
2 


*• 

|; P»a 


»P» 1 


»l 


( 34 ) 


Using these approximations, (25) becomes 


Pr{n>77l^} = | 




f^ + /^ exp L 

1 2* V 2* J l 


(Vfr -Va) 2 j 


(35) 


or, from (26), 


Pr{f) > Tjld^} = 


N-y j 

^ 2 -|5f 

N + yj 

^V 2 -|5| 2 


V' 4 , 
+ 



^ 2 -|5f 

N-^l 

^N 2 - 15| 2 


\ 1/4> \ 


> ) 




2An- 


N n 


+ lexpj — —(N - |5|)| 

!|5I(A^-|<5|) I y \ 2 N, 1 '7 


(36) 


For N = 3 and M = 2, 151 =1 from (32). Thus, (33) becomes 





(37) 


Comparing (37) with (31), we observe that the factor in front of the term in 
brackets in (37) represents a bound on the improvement in performance 
obtained by increasing the memory of the decision by one symbol interval 
from N = 2 to N = 3. 


3.4 General Asymptotic Results 

In the general case for arbitrary N, the dominant terms in the bit error 
probability occur for the sequences that result in the minimum value of N 
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151. One can easily show that this minimum value will certainly occur for 
the error sequence A0 having N-l elements equal to the correct sequence A0 
and one element with the smallest error. Thus, 


min (N-l 51) = N-\N- 1 + e ***-**>nu» 

1 1 

= N-V(N- l) ! + (IV-lX2-<&,) + l 

= N-I«U. 

(38) 

where 

d 2 . - 4sin 2 ^^ — QQjshm = 4 sin 2 — (39) 

rain 2 M 

Also note that for 151 = 151 max , (26) reduces to 



= -^~ N±2VN-lsin— 

2N 0 L M J 


(40) 


Thus, the average bit error probability is approximately upper bounded by 




1 

(N-l)log 2 A/ 

r 


^(M) 




1 



PLTO 


ex pi 


2 N, 


;(w-|5L)| 


(41) 


where w(u,u) corresponds only to those error sequences that result in 1 5 1 max- 


For the binary case (M = 2), we have from (39) that d ^ = 4 and hence 
N-1 5l nux = 2. Similarly, it is straightforward to show (see Appendix B) that the 
sum of Hamming distances required in (41) is given by 



12 


£vv(w,u) = 


[U N = 2 


Thus, (41) simplifies to 


N> 2 


(42) 


n<- 


2 f /OYI 4-ill 

rri,Viv-2j[2 p i n 0 jj 

1 


which is the generalization of (37) for arbitrary N 2 . 


(43) 


Eq. (43) has an interesting interpretation as N gets large. Taking the limit 
of (43) as N -* oo, we get 



which can be expressed in terms of the asymptotic expansion of the 
complementary error function. 



(45) 

(46) 


For coherent detection of binary PSK (BPSK) with differential encoding and 
decoding, the bit error probability performance is given by [2: Chap. 5] 



which has an asymptotic upper bound identical to (46). Thus, as one might 
expect, the performance of multiple symbol differentially detected BPSK 
approaches that of ideal coherent detection BPSK with differential encoding 
in the limit as the observation interval (decision memory) approaches 

infinity. 


2 Note that (43) is not valid for N = 2 since in that case |8 1 ^ * 0 Isee (38) and (39)1 and 
thus the inequalities in (34) are not satisfied. 
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A similar limiting behavior as the above may be observed for other values 
of M. In particular, it can be shown (see Appendix B) that for M > 2 and a 
Gray code [2] mapping of bits to symbols, the sum of Hamming distances 
corresponding to 1 51 max is given by 


'£w(u,u) = 


[4(iV-l); N> 2 
1 2; N = 2 


Using (48) in (41), we get (for N > 2) 

2 




OoMj 


N + SI 

Hn» 


|«L<Af-|«L>. 


where, from (38) and (39), 




I SL.= J(N - 1) 2 + 2 (N - 1X1 - 2 Sin 2 -^;) + 1 

M 


For N = 2, the upper bound on bit error probability becomes 




1 


(\og 2 M)J2njf 


COS 

n 

2M 

n 

K 

sin 

\ cos-— 

2 M 

\ M 


xexpj- 


2E, . 2 n 

— t sin 

N 0 2M 


As N gets large, -> N - 2sin 2 -^ and (49) reduces to 


P„< 


1 


(lo g! M)J)t#-sin ! ^ 


I E , • 2 n 

ex P i 77 sm 77 


1 


N 0 M I log 2 M 


-erfc 


(48) 


(49) 


(50) 


(51) 


'fir. 

— sin — 
Mj 


N n M 


( 52 ) 
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which is identical to the asymptotic bit error probability for coherent 
detection of MPSK with differential encoding and decoding (see [2: Eqs. (5-91), 
(5-92) and (5-113)1)3. 


For example, for QPSK (M=4), the symbol error probability is given by 
[2: Eq. (5-115)] 


P. = 2erfc 


2N. 




-2erfc i 


oy 


J 2N 0J 


+ erfcl 


pK) 4 [pN 0 






(53) 


Since for a Gray code bit to symbol mapping 




(54) 


log 2 M 2 

then (54) together with (53) has an asymptotic upper bound identical to (52). 


Figures 4, 5, and 6 are illustrations of the upper bounds of (49) and (51) for 
M = 2, 4, and 8, respectively. In each figure, the length (in MPSK symbols) of 
the observation interval, N, is a parameter varying from N = 2 (conventional 
MDPSK) to N = oo (ideal coherent detection). Also indicated on the figures 
are computer simulation results corresponding to the exact performance. We 
observe that, for example, for binary DPSK, extending the observation 
interval from N = 2 to N = 3 recovers more than half of the E^/Nq loss of 
differential detection versus coherent detection with differential encoding. 
For M = 4, the improvement in E^/Nq performance of N = 3 relative to N = 2 
is more than 1 dB which is slightly less than half of the total difference 
between differential detection and coherent detection with differential 
encoding. 


4.0 Application to Trellis Coded MDPSK 


In this part of the report, we extend the idea of multiple differential 
detection of MPSK to trellis coded modulations (TCM). We will show that a 


3 It should be noted that the result in (52) can be obtained by observing that, for large N, (40) 
satisfies Vb » Vb - Va > 0. In this case, (25) can be approximated by [5: Appendix A] 

Pr{r/ > r/ld£} = Ul - Q(s[b , Vfl ) + fi(Va, Vft)] s ^erfcf ^ ^ ~ ^ 




V2 


Using this relation in (22) gives the asymptotic bit error probability in (52). 
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combination of a multiple trellis coded modulation (MTCM) [7] with 
multiplicity equal to N-l combined with multiple 4 symbol differential 
detection can potentially yield a significant improvement in performance, 
even for small N, over that corresponding to conventional trellis coded 
MDPSK. 

The analysis technique that will be employed to obtain upper bounds on 
the bit error probability performance of the system is equivalent to that used 
in [8] to assess the performance of conventional trellis coded MDPSK on a 
fading mobile satellite channel. In fact, it will be shown that the results 
obtained here have an interesting similarity to those obtained in [8] once an 
association is made between the squared Euclidean distance measure per 
trellis branch for conventional differential detection and the equivalent 
distance measure per trellis branch for multiple symbol differential detection. 

5.0 System Model 

Figure 7 is a simplified block diagram of the system under investigation. 
Input bits occurring at a rate are passed through a rate nk/(n+l)k multiple 

trellis encoder (k is the multiplicity of the code) producing an encoded bit 
stream at a rate R s = [(n+l)k/nk]R^. Next, the encoded bits are divided into k 

groups of n+1 bits each and each group is mapped into a symbol selected 

from an M = 2 n+ ^ - level PSK signal set according to a set partitioning 
method for multiple trellis codes [7] analogous to that proposed by 
Ungerboeck [9] for conventional (unit multiplicity) codes. Since the MDPSK 
symbol rate is R^/n, it is reasonable, from a conservation of bandwidth 

standpoint, to compare the performance of this system to an uncoded M = 2 n 
level DPSK system with the identical input bit rate. 

At the receiver, the noise-corrupted signal is differentially detected and the 
resulting symbols are then inputted to the trellis decoder which is 
implemented as a Viterbi algorithm. In selecting a decoding metric, a tradeoff 
exists between simplicity of implementation and the optimality associated 
with the degree to which the metric matches the differential detector output 
statistics. 

4 One must be careful not to confuse the multiplicity of the trellis code with the 
multiplicity of the differential detection scheme although, as we shall soon see, the two 
are indeed related. 
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For the case of uncoded MDPSK, a metric based on minimizing the 
distance between the received and transmitted signal vectors is optimum in 
the sense of a minimum probability of error test. The specific forms of this 
metric for conventional and multiple differential detection were described in 
Section 2.0. For conventional trellis-coded MDPSK, the metric takes on the 
form of a minimum squared Euclidean distance metric. For multiple symbol 
detection of MTCM, the form of the metric is quite different . Nevertheless, as 
we shall soon see, by a suitable modification of the multiple trellis code 
design, the appropriate metric can be converted once again into a minimum 
squared Euclidean distance metric. The so-called "equivalent" multiple 
trellis code that results from this modification then becomes the key tool used 
for analyzing the performance of the system. 

6.0 Analysis Model 

We denote a coded symbol sequence of length N s by 

x = {x v x 2 ,---,x N ) (55) 

where the kth element of x, namely, x^, represents the transmitted MPSK 
symbol in the kth transmission interval and, in general, is a nonlinear 
function of the state of the encoder and the nk information bits at its input. 
Before transmission over the channel, the sequence x is differentially encoded 
producing the sequence £. In phasor notation, sfc and s^+i can be written as 

s k = y(2Pe j * k 

s M = s k x M = V27V (# * + ** +,) = V2 Pe» M (56) 

where E s = rE^ is the energy per MDPSK symbol and 

x k = e iMk (57) 

is the phasor representation of the MPSK symbol A0* assigned by the mapper 
in the kth transmission interval. 

The corresponding received signal in the kth transmission interval is 
given by (2) with s^ now defined as in (56) and the noise sample, n^, defined 
as before with variance given by (3). 
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Consider now a received sequence of length N s and assume that 0^ is 
independent of k over the length of this sequence, i.e., 0^ = 0. Then, the 

received sequence r is expressed as in (4) where r^, sfc, and n^ are, 
respectively, the kth components of the N s -length sequences r, s, and n. 

Since the detection scheme will be independent of 0, we can furthermore set 
0 = 0 without any loss of generality. 

To apply the notion of multiple symbol differential detection to trellis 
coded MPSK, the decision statistic of (15) must be associated with a branch in 
the trellis diagram. To do this, we construct a multiple trellis code of 
multiplicity k = N-l . (The procedure for designing this code will be discussed 
later on in the publication.) Thus, we can envision the transmitted sequence, 
x, of (55) as being partitioned into B = N s /k = N S /(N-1) subsequences 5 , i.e.. 


x = (x l '\x l2 \....,x' B) ) 


(B)\ 


(58) 


with each subsequence xw = (xjj, xj2, — /*ik) representing an assignment to a 
trellis branch. Similarly, a received sequence, r, of length N s is associated 
with a path of length B branches in the trellis diagram. Once this 
association is made, computation of bit error probability for the system 
follows along the lines of the approach taken in [7]. The details of the analysis 
are presented in the following sections. 

7.0 Derivation of Pairwise Error Probability Bound 

To find an upper bound on the average bit error probability performance of the 
system, we must first find the pairwise error probability which represents the 
probability of choosing the coded sequence x = (x l ,x 2 ,...,x Nj ) = (x°\ x (2 \...,x (B) ) instead 
of x = (x l ,x 2 ,...,x Nt ) = (x m ,x {2) ,...,x iB) ). Letting qj denote the maximum-likelihood 
metric for the correct data phase sequence on the ith trellis branch and computed in 
accordance with (15), then the pairwise error probability is given by 

A<f>\ (59) 




1*1 


»*1 


5 Since N s is arbitrary, we can choose it such that N S /(N-1) is integer. 
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Here r), denotes the metric computed for the data phase sequence associated with 
the ith trellis branch of the incorrect path, A0 denotes the transmitted data 
phase sequence for the correct path, and B is the length (in branches) of the 
correct and incorrect paths. Also, since the channel is additive white Gaussian 
noise (AWGN), the summations in (59) represent sums of independent random 
variables. 

An exact evaluation of (59) in closed form is difficult if not impossible. At 
first glance it might appear that the decision variable in (59) is a special case of the 
quadratic form given in [10: Eq. (4B.1)] and thus the pairwise error probability 
would be given by [10: Eq. (4B.21)]. Unfortunately, however, the development in 
[10: Appendix 4B] requires that the second central moments of the complex- 
valued random variables r}., ff ( ; i = 1,2,.. .,B be independent of i. While indeed 
the variances of these random variables are independent of i, the covariance of 
TJi and i) i is directly proportional to 8j (see Eq. (A-5)) defined by 


N-i «-2 


N-n-2 


N -i > N-l > X 

" =0 = X e " =0 


*=0 


*=0 


(60) 


and thus depends on i. In (60), the data phases that appear in the exponent are 
the elements of A0 (l) and A0 (,) which denote the ith subsequences of A0 and Atp_ , 


respectively. 


In principle, then, there are two approaches one can take to evaluating ( 59 ). 
The first is to derive a Chemoff bound [11] on (59) using a method similar to that 
taken in [8]. This method has the advantage (as we shall soon see) of enabling 
the upper bound on bit error probability to be obtained using the transfer 
function bound approach applied to trellis coded modulations in [7]. It has the 
disadvantage (typical of Chemoff bounds) of resulting in a loose upper bound on 
bit error probability. 


The second approach is to try to directly approximate (rather than upper 
bound) the result in (59) for large signal-to-noise ratio (SNR). In Section 3.0 we 
saw for the uncoded that, despite the fact that we were able to exactly evaluate the 
pairwise error probability (equivalent to evaluating (59) for a path consisting of a 
single trellis branch, i.e., B = 1), we eventually approximated that result for large 
SNR and showed that the results agreed extremely well with those obtained 
from simulation. Also, the asymptotic (large SNR) form of the bit error 
probability expression had the advantage of allowing direct comparison with an 


19 


analogous expression for conventional (no multiple symbol observation) 
MDPSK. The disadvantage of this approach is that the expression for the 
approximate pairwise error probability that results cannot be put in a form that 
allows evaluation of the upper bound on bit error probability by the transfer 
function bound approach. Thus, in this part of the report, we shall consider 
both approaches because of their respective merits. 

7.1 A Chemoff Bound on Pairwise Error Probability 

A Chemoff upper bound [6] on the pairwise error probability of (59) is 
evaluated in Appendix C with the result 


P(x ->x)<, 






£ A(l-^)[iV 2 -|5 i | 2 ]) 

cxp< 

N 0 1- A 2 

* 2 -|5. 2 

J 


l-A 2 [iV 2 - 

-Itf] 



where 5i is given by (60). 


(61) 


The expression in (61) bears a striking resemblance to that which characterizes 
the pairwise error probability of conventional differential detection of trellis 
coded MDPSK. In particular, letting p n = 1 (no fading) and X, = 2Xq in [8: Eq. (25)], 
we get 


a ex P 

P(x i) < fl 

7=1 


E, A(1-2A fe-jc,| 2 
N 0 

i-W-4 


(62) 


Note that for N = 2 (conventional differential detection) and M = 2 (binary 
DPSK), 


<S, = 0 (63) 

in which case, (61) agrees with (62). 


Comparing (61) with (62) we observe that, for multiple symbol differential 
detection, the equivalent squared Euclidean distance measure per trellis branch 
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is - 1 5|2. We now examine in detail the algebraic structure of - 1 8|2. In 
particular, we shall show how - 1 8|2 can be expressed as a squared Euclidean 
distance of an "equivalent" trellis code with larger multiplicity. Once this 
association is made, this "equivalent" trellis code then serves as the 
mathematical tool by which the foregoing analysis in [3] can be directly applied 
here. 

7.1.1 The Construction of an Equivalent Multiple Trellis Code with Squared 
Euclidean Distance Measure 


Here we show how to take a trellis code of multiplicity N-l and distance 
measure - 1 8|2 and construct an equivalent (in performance) code with 
larger multiplicity but a squared Euclidean distance measure. For simplicity 
of explanation, we will start with the case N = 3 which corresponds to only 
one additional symbol of observation relative to conventional differential 
detection. 


From (57) and (60), we have that 
N 2 - 1<5,| 2 = N 2 - 


= N 2 -|l + 

For N = 2, (64) simplifies to 


N - 1 N-h-2 I 2 

s iwusd 

«*0 I 

N-2 N-n - 2 I 2 

£ n*jed 

*=0 *i *=0 I 


(64) 


N 2 - |5j 2 — 4 — jl -+- ) *| 2 (65) 

Since there is only one MPSK symbol per trellis branch, i.e., multiplicity equal 
to one, we can simplify the notation in (65) to 


N 2 - \8\ 2 = 4 - |l + = 4 - (l + + 2Re{x i i;}) 

= 2 - 2 Re { } = |^, - ( 66) 

Thus, for conventional differential detection, - 1 8|2 yields the squared 
Euclidean distance measure as one would expect. 


For N = 3, (64) becomes 


21 


« 2 - fcf = 9 - 11 + 

= 9 - (3 + 2 Re{^>x<‘> + 2Re{x“ jt*2T} + 2 Re{x<‘>if}) 

= 6- (2 - |*M, - «/) - (2 - |xf2, - i“f) - (2 - |*» - iff) 

= M?, - iaf + |4° - *?f + l*M- - «-.f 

(67) 

The first two terms in (67) represent the squared Euclidean distances 
associated with the two symbols assigned to the if/i trellis branch and as such 
their sum would be the squared Euclidean distance for this branch. The third 
term in (67) can be interpreted as follows. Note that the product of two MPSK 
symbols, say and x m , is indeed another MPSK symbol whose phase is the 
modulo M sum of the phases of x^ and x m . Thus, if the MPSK symbols are 

represented by their equivalent M-ary numbers, e.g., el^nt/M m/ then the 
third term in (17) represents the squared Euclidean distance associated with a 
symbol which is the modulo M sum of the first two symbols assigned to that 
trellis branch. 

As a simple example of the above construction, consider the multiple 
trellis code illustrated in Figure 8 which has multiplicity k = N - 1 =2. This 
code is the optimum 2 state, k = 2, rate 2/4 trellis coded QPSK designed for the 
AWGN. Computing the performance of this code with a three symbol 
multiple differential detection scheme, the distance measure - 1 8|^ is then, 
according to the above, mathematically the same as computing the 
performance of the k = 3 trellis code in Figure 9 and a squared Euclidean 
distance measure. Note that the third symbol assigned to each branch in 
Figure 9 is the modulo 4 sum of the first two symbols on each path of that 
branch. 

The above technique for generating the equivalent code can be generalized 
to arbitrary N. In particular, from (64) it is straightforward to show that 6 


6 For simplicity of notation, we omit the "i's" on the variables with the understanding 
that we are dealing with the i th branch. 
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N- 2 


N - 3 JV-2 


W]=0 «i=0 j«2-i 

tf-4 N - 3 AT-2 2 

X | * -"l”** - *2^k - "3 *3| 


iq=0 «2= 1 *3=2 

+ • "kt**-r • ' x k-s + 2 ~ *t*k- r • •-^t-jv+J 2 


( 68 ) 


Once again the first term in (68) represents the squared Euclidean distance 
associated with a branch in the N-l multiplicity trellis code whereas the 
remaining terms represent the squared Euclidean distances of the additional 
symbols that must be assigned to each branch in accordance with modulo M 
sums of the previous symbols. 

It is important to emphasize that the trellis code of Figure 9 (or, more 
generally, the one that would be constructed from (68)) is strictly a 
mathematical tool that is useful for performance analysis. The actual trellis 
code that would be used for an N = 3 multiple differential detection of rate 2/4 
trellis coded QPSK would be a multiplicity 2 code such as that illustrated in 
Figure 8. 

7.2 An Asymptotic (Large SNR) Evaluation of Pairwise Error Probability 

As previously mentioned, finding an exact closed form expression for the 
pairwise probability of (59) is difficult due to the dependence of the covariance 
of tj ; and f), on the summation index i. Nevertheless, if one is willing to 
settle for asymptotic (large E s /No) results, then it is possible to define a 
procedure which will allow approximate evaluation of pairwise error 
probability. 


To see how this comes about we first recognize that the decision variable is 
a special case of the quadratic form discussed in [10: Appendix 4B]. Thus, 
following the approach taken there, it is straightforward to show that for the 
specific case at hand here where the second central moments depend on i, 

[10: Eq. (4B.9)] generalizes to 
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where 

« 2 , = (2P)[N , -|5,f] 


V< = 



(70) 


with 5j defined in (60) and e > 0 a parameter that can be selected, as we shall 
soon see, for convenience in defining the contour path of integration. At 
this point, we cannot proceed further with the procedure taken in 
[10: Appendix 4B] due to the dependence of vj on i. 


A technique for evaluating (70) in the limit of large PT/Nq = E s /No is 
discussed in Appendix D and is based on a variation of the method of 
stationary phase for analytic functions [13]. In order to apply this technique to 
(70), we first rewrite it in the normalized form 


P , x _> i) L r"' irr_J_ ex jA.y JMmzIK&L, 


(71) 


where 


A-w 

if 



(72) 


and e'= e/v ; . where i m in * s the va ^ ue of i that results in the minimum 
hnin 

value of 1 I , i.e., the trellis branch along the error event path with the largest 
equivalent squared Euclidean distance N2 - 1 1 2. The integral in (71) is now 
in the form of (D-l) where 


’ h «,v+i 


(73) 

Before considering the general result for an error event path with an 


24 


arbitrary number of branches B, we shall show that, for B = 1 (a one branch 
path), this approach gives the identical result to (36) for the uncoded case. In 
particular, for B = 1, we have from (72) that = 1, and thus letting 5 = 8|, the 

appropriate 7 critical point (value of z where f (z) = 0) of (73) is given by 


0 ’yN + |5| 1 


(74) 


Also, note that since the integration contour in (71) is selected to pass through 
the critical point of (74), then 


e ._ 


(75) 


The real and imaginary parts of f(z) evaluated at the critical point of (74) 
are given by 


Also, 


0 1 (z o ) = Re{/(z o )} = O 
0 2 (z o ) = Im{/(z 0 )} = 


/’(* o) = 


y (N+l£j) 2 

2 \5\ 


gU 0 ) = -j 


W+]g| N+j6 

2 I 5 I 


(76) 


(77) 

Finally, substituting (76) and (77) into (D-13) gives the desired result 


P(x -> i) = Prjrj > ri\A<p} = 



' N + \S\ " 

\SKN-\8\)j 




which agrees with (36). 



(78) 


7 Another critical point of f(z) occurs at 


• >+| 5 | 

N-\d\ 


Zo -J 


>j 


However, this point results in the integration contour passing through one of the 
singularites of g(z) at z = j and thus, in accordance with Appendix D, is not allowable. 
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For an error event with an arbitrary number of branches, B, we reason as 
follows. Since the parameters {^; i = in (72) have been defined such 

that their values are all less than or equal to unity, then from (73) we observe 
that, aside from the singularity of g(z) at z = 0, all of the singularities of f(z) 
and g(z) lie on the imaginary z axis in the interval j < z < °°. Furthermore, it 
is shown in Appendix E that there exists a unique critical point of f(z) in the 
interval 0 <, z < j. Denoting this critical point by zq = jyg (unfortunately, one 
cannot, in general, find a closed form expression for zq analogous to (74) and 
thus one must resort to numerical evaluation), then the integration contour 
in (71) should again be selected to pass through this point, i.e., choose e' = yo- 
Furthermore, the contour should be tailored so that it does not pass through 
any of the other critical points of f(z) should they occur at values of z = a + je'; 
a # 0. 

From (73), we observe that for zq = jyo, f(zg) and g(zg) are purely 
imaginary, i.e.. 


g(m) = -J 


SU 

If i-J— 

Ah-$k 2 


i=l * hiJ 


=-70 2 (> b) 


(79) 


Also, the second derivative of f(z) evaluated at z = zq is also purely imaginary 
and given by 




1=1 


M^ 2 ) 3 


=j<t> a(3b) 


(80) 


Then, using the results of Appendix D, in particular, Eq. (D-13), the pair- 
wise error probability of (71) is asymptotically approximated by 


P(x ->x) = 


r 0 2 (3b) ^ 


2 n- 




e x P --to) 


(81) 


N n 


which is of the same functional form as (78). 
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8.0 Evaluation of an Upper Bound on Bit Error Probability 

An upper union bound on the average bit error probability is obtained 
from the pairwise error probability as 

P b ^ ^^a(x,x)p(x)P(x -> x) (82a) 

l.xec 

where a(x,|)) is the number of bit errors that occur when x is transmitted and 
x is chosen by the decoder, p(x) is the a priori probability of transmitting x and 
C is the set of all coded sequences. If the pairwise error probability is upper 
Chernoff bounded as in Section 7.1, then one must in addition optimize over 
the Chernoff parameter, in which case (82a) becomes 

P b <, rnin ^^a(x,x)p(x)P(x -> x) (82b) 

xMc 

8.1 Evaluation bv the Transfer Function Bound Approach 

When the pairwise error probability is upper Chernoff bounded, then an 
efficient procedure for evaluating (82b) is the transfer function found 
approach applied to multiple trellis coded modulations in [7]. In particular, 
the trellis codes are represented by a pair-state transition diagram [14]. Each 
pair-state (s if s t ) corresponds to a pair of states s ( and in the trellis diagram. 
Thus, a transition between pair-states ($,,£) and (s M ,s M ) in the transition 
diagram corresponds to a pair of transitions in the trellis diagram, i.e., s i to s M 
and ^ to s i+1 . Associated with each of these trellis diagram transitions are the 

k' MPSK symbols 8 corresponding to a sequence of nk input bits (an 
information symbol) to the multiple trellis encoder. Thus, the transition 
between two pair-states in the transition diagram is characterized by a 

function of the squared Euclidean distance |z,- - z,j 2 ; i = 1,2,.. .,k' between the 
corresponding k' MPSK output symbols^ and the Hamming distance £2 
between the corresponding input bit sequences. 


®Here, k' refers to the extended multiplicity of the equivalent code as discussed in the 
previous section. 

®Here z ( - and z { denote the correct and incorrect MPSK symbols assigned to a trellis 

branch in the equivalent code with extended multiplicity. As such, the first k of these 
symbols will be identical to the x/s and Jt.’s, respectively, whereas the remaining 

symbols represent modulo M sums of these first k in accordance with the construction 
procedure discussed in Section 8.1. 
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Based on the above discussion, each branch between pair states in the 
transition diagram has a gain G of the form 

al 


exp 




Ad-AAoXfc-zj 


k' 

1 

i= 1 


N n 


1 - - 4 


;=i 






(83) 


i=l 


where I is an index, nk is the number of bits input to the trellis encoder in 
each transmission interval, and 


/(*, A)= 


4A(1-W)» 

& o i-* 2 * 


l-X*x 


(84) 


with D the Bhattacharyya parameter [11,14] defined by 


_ 4 

D = exp 


4AT 


(85) 


o J 


Also in (83), the summation accounts for the possibility of parallel paths 
between states in the trellis diagram. The transfer function (the sum of all 
possible path gains) of the transition diagram is denoted by T(D,I) and, by 
comparison with (82b), the upper bound on average bit error probability is 
given by 


P„ Z min— — T{D,I)\. . (86) 

* i nkdl l/=1 

8.1.1 A Simple Example 

To illustrate the foregoing theory, we consider, as a baseline, the simple 
case of a 2 state, rate 1/2 trellis coded QPSK designed for the AWGN. The 
appropriate trellis diagram for such an Ungerboeck design is illustrated in 
Figure 10. If indeed we are to detect this coded modulation differentially with 
N symbol observation, then we must first construct a multiple trellis code 
with multiplicity k = N-l which would have identical performance as the 
original k = 1 code of Figure 10 if conventional differential detection (N = 2) 
were used for both. 10 

10 Note that the construction of such an equivalent (same performance) multiple trellis 
code from a unit multiplicity (Ungerboeck-type) trellis code can always be performed 
whereas the reverse is not necessarily true. Also note that the multiple trellis code does 
not require a different encoder than that used for the unit multiplicity code. 
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Hie particular case we shall consider here is N = 3. Thus, we need to 
transform the trellis diagram of Figure 4 into an equivalent multiple trellis 
diagram with multiplicity k = 2. The procedure for doing this is to perform a 
particular type of Cartesian product of the transition matrix of the code in 
Figure 10 with itself. Letting 



represent the transition matrix of a 2 state code with no parallel paths, i.e., xjj 
is the MPSK symbol assigned to the branch corresponding to the transition 
from state "i" to state "j", then the appropriate transition matrix for the 
equivalent multiple trellis code is 



From (75), we see that each branch in the equivalent trellis code is 
characterized by two parallel paths with code symbol assignments in 
accordance with the 2-tuples in T^). 


Applying (87) to Figure 10, we get 


T = 


'0 2 ' 
1 3 


( 89 ) 


and thus, from (75), the transition matrix of the equivalent code is given by 



which generates the trellis diagram illustrated in Figure 11. Note that this 
trellis is different than that illustrated in Figure 8 for the same rate 2/4 coded 
QPSK. As such. Figure 11 does not represent the optimum multiplicity two 
code. 


Figure 12 illustrates the pair-state transition diagram corresponding to 
Figure 11 and has the transfer function 
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T(D,I) = 


Aac 


a- 


b = 


c = 


exp 

I 

^2b 

£ k 4A(l-2A)"| 

N 0 1-4A 2 J 

2 

exp 

/ 

1-4A 2 

E b 2A(1-2A)^ 

N 0 1-2A 2 J 

2 

exp 

1 

1-2A 2 

f E b 2A(1-2A)"| 

l 1-2A 2 J 


1-2A 2 


b 

I 


(91) 


where f(x,A) is defined in (84) with N = 2. Differentiating (91) in accordance 
with (86) gives the desired expression for the upper bound on bit error 
probability, namely. 


D „ . /(4,A)/(2,A) . , 

P. < min^ — — ~r- = min 

* (1-/(2, A)) 2 * 


^l-2A n 

1-4A 2 


exp 


E b 2A(3-8A 2 ) 
l Nq (1+2A)(1-2A 2 ) 


1 - 2A 2 - exp - 


E k 2A (1 - 2A ) 
N 0 (1-2 A 2 ) 


YV 


JJ 


(92) 


The upper bound on of (92) is plotted in Figure 13 versus E^/Nq. 


To apply multiple symbol detection to the trellis of Figure 11, we must first 
convert it, in accordance with Section 8.1, to an equivalent trellis diagram 
with a Euclidean distance measure for all its symbols. This mathematically 
equivalent trellis diagram is illustrated in Figure 14 for N = 3. Figure 15 
illustrates the pair-state transition diagram corresponding to Figure 14 and 
has the transfer function 


r(D , /)=2c+ _±|_ 

a = Ul + / 2 )/(8,A); b = j[(2 + /)/( 4, A ) + //(8, A )] 

2 4 

c = |//(8,A); <f = i[(2/ + 7 2 )/(U) + / ! /(4,*)] 

where f(x,A) is defined in (84) with N = 3. Differentiating (93) in accordance 
with (86) gives the upper bound on bit error probability as 
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P^mmi/(8,A)| 


1 + 


1 1/(4, A ) + 5/(8, A) - 5/(4, A )/(8, A ) + ^-/ 2 (4, A ) - |/ 2 (8, A ) 


2 1-|/(8,A)-|/(4,A) 


(94) 


The upper bound of (94) is superimposed on the results from (92) in Figure 
13. Also included in this figure are simulation results corresponding to the 
exact performance of the system for N - 2 and N = 3. 


We observe from Figure 13 that the curves corresponding to the Chemoff 
bounds of (92) and (94) are quite close to one another. The reason for this is 
that, for N = 3, the bound is much looser than it is for N = 2. To understand 
the reason for this, we compare the exact result (or its asymptotic 
approximation) obtained for uncoded multiple symb ol DPSK with what 
would have been obtained by using a Chernoff bound on the pairwise error 
probability. In particular, if instead of the asymptotic approximation of (36) 
the pairwise error probability were Chemoff bounded (as has been done here 
in the coded case), then one would obtain a relation analogous to (61) but 
without the product since, for the uncoded case, the number of branches in an 
error event, B, is equal to one. If then one were to minimize only the 
exponential term in (61) with respect to X (global minimization is difficult to 
accomplish in closed form), then the optimum Chemoff parameter becomes 
Xopt = 1/(N + 181) which when substituted in the Chemoff bound gives 

(95) 

Comparing (95) with (36), we observe that both the approximation of the exact 
result and the Chemoff bound yield the same exponent; however, the 
Chemoff bound does not produce the inverse square root of symbol energy- 
to-noise ratio behavior which is all important in distinguishing the 
performance of conventional from multiple symbol differential detection. 
Thus, using the Chemoff bound of (95) rather than the approximation of the 
exact result as in (36) to calculate the union bound on bit error probability for 
multiple symbol differential detection of uncoded MPSK would also result in 
a loose upper bound as we have observed in the coded case. 
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For equivalent trellis diagrams with parallel paths (as is the case in the 
current example - see Figure 14), we can improve upon the above as follows. 
Consider the terms in (82b) due to the parallel paths. Since parallel paths are 
one branch in length, then for these terms we can use the exact (or large SNR 
approximate) result for P(x — » x) from the uncoded results in Section 3. In 

terms of the evaluation of by the transfer function method, what we do in 
effect is to subtract the portion of the transfer function due to the parallel 
paths and apply the Chernoff bound to only the remaining portion of the 
transfer function. Mathematically speaking, this is equivalent to rewriting 
(86) as 

/>S/V. + mn^|-[r(A/)-7i(D,/)] J , | (96) 

where P^q is the portion of P^ contributed by the parallel paths; this portion is 
evaluated by the exact (or large SNR approximate) result, and Tq(D,I) is the 
part of the transfer function due to these parallel paths. The amount of 
improvement obtained using (96) instead of (86) will be significant when the 
terms due to the parallel paths dominate the error probability, i.e., they yield 
the minimum equivalent squared Euclidean distance (N^ - 1 5 1 2) over all 
error event paths. 

For the example under consideration, the term "2c" in (93) results from 
the parallel path in Figure 14, i.e., Tq(D,I) = 2c. This, in turn, contributes the 
term (l/2)f(8,X) in (94) which corresponds to the minimum equivalent 
squared Euclidean distance N^ - 181^ = 8. Subtracting this term from (94) and 
replacing it by the asymptotic (large SNR) approximation [see (36)] 



where for this example, N = 3 and 



then, from (94), (96), and (98), we get the approximate tighter upper bound 
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1 1/(4, A ) + 5/(8, A ) - 5/(4, A )/(8, A ) + |/ 2 (4, A ) - 1 / 2 ( 8, A ) 
2[l-|/(8,A)-i/(4,A)] 

(99) 

The bound of (99) is superimposed on the results in Figure 13 and is seen to 
yield an improvement over that corresponding to (94). 

8.2 Evaluation Using the Asymptotic Approximation to Pairwise Error 

Probability 

The upper bound of (82a) together with (81) has been evaluated for the 
example under consideration. In performing this evaluation, we have kept 
only those error event paths which yield a signif icant contribution to the 
sum. The results are also illustrated in Figure 13 and are seen to represent an 
improvement over the upper bounds determined in Section 8.1. However, 
since for this example, the trellis contains a parallel path which will dominate 
the error probability performance at high SNR, then the upper bound of ( 99 ) 
is quite close to that obtained here using the asymptotic approximation to 
pairwise error probability. For trellises that do not contain parallel paths, the 
approach leading up to Eq. (99) cannot be used and thus one must employ 
either the loose Chernoff bound or the asymptotic approximation method. 

The same upper bound of (82a) together with (81) can also be used for 
conventional (N = 2) differential detection of rate 1/2 trellis coded QPSK. For 
the example in Section 8.1.1, the result is illustrated in Figure 13 and is again 
seen to agree quite well with the comparable simulation results. 

8.3 Another Example 

Figure 16 illustrates error probability performance results obtained by 
simulation for a 16 state, rate 2/3 trellis coded 8PSK using conventional 
(N = 2) and multiple (N = 3) symbol differential detection. This code, which 
is optimum on the AWGN, has the transition matrix [9: Fig. 7] 


+ min^/(8,A)| 
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T = 


0 4 2 6 


15 3 7 


2 6 0 4 


4 0 6 2 


6 2 4 0 


4 0 6 2 


3 7 15 


5 17 3 


7 3 5 1 


5 17 3 


6 2 4 0 


7 3 5 1 


0 4 2 6 


15 3 7 


2 6 0 4 


( 100 ) 


3 7 1 5_ 

We observe from this figure that in going from N = 2 to N = 3, an 
improvement in E^/Nq performance of about 0.75 dB is obtained for the 
range of bit error probabilities illustrated. 


9.0 Conclusions 

We have demonstrated a multiple symbol differential detection technique 
for MPSK which is based on maximum-likelihood sequence estimation 
(MLSE) of the transmitted phases rather than symbol-by-symbol detection. 

The performance of this multiple symbol scheme fills the gap between ideal 
coherent and differentially coherent detection of MPSK. The amount of 
improvement gained over differentially coherent detection (two symbol 
observation interval per decision) depends on the number of additional 
symbol intervals added to the observation. In the limit as the observation 
interval approaches infinity, the performance approaches that of ideal 
coherent detection with differential encoding. Practically, this limiting 
performance is approached with observation times (decision memory) only 
on the order of a few additional symbol intervals. Thus, even in situations 
(e.g., benign environments) where one would ordinarily not turn to 
differential detection, it might now be desirable to employ multiple symbol 
differential detection for reasons related to simplicity of implementation. For 
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example, the acquisition and maintenance of a locked carrier tracking loop as 
required in a coherent detection system is not needed here. 

As for the uncoded case, the use of multiple symbol differential detection of 
trellis coded MPSK can also offer an improvement in error probability 
performance over conventional (two-symbol observation) differential 
detection of the same coded modulation. Again only a slight increase in the 
length of the observation interval is necessary to demonstrate a significant 
improvement. 
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Figure 5. Bit Error Probability Versus E^/Nq for Multiple Differential Detection 
of MPSK; M=4 
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Figure 6. Bit Error Probability Versus E^/Nq for Multiple 
Differential Detection of MPSK; M=8 
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Figure 7. Block Diagram of the Trellis Coded MDPSK System 
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Figure 9. Trellis With Squared Euclidean Distance 

Measure That is Mathematically Equivalent to 
Figure 8 



Figure 10. Trellis Diagram for Conventional Rate 1/2 
Trellis Coded QPSK 
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Figure 11. Trellis Diagram for Rate 2/4 Multiple Ti 
Coded QPSK Equivalent to Figure 10 



igram for Trellis Diagram of Figure 11 


46 



Figure 13. Analytical and Simulation Results for Bit Error Probability of 2-State, 
Rate 1 /2 Trellis Coded QPSK With Conventional (N=2) and Multiple 
(N=3) Symbol Differential Detection 
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0 , 0,0 0 , 0,0 



Figure 14. Trellis With Squared Euclidean Distance 

Measure That is Mathematically Equivalent to 
Figure 11 
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[I] 


16 State TCM 
Interleaving Size = 512 symbols 
Interleaving Depth = 1 6 branches 


E b /N 0 ,dB 


Figure 16. Simulation Results for Bit Error Probability of 16 State, Rate 2/3 
Trellis Coded 8PSK With Conventional (N=2) and Multiple (N=3) 
Symbol Differential Detection 
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Appendix A 

Evaluation of the Pairwise Error Probability 


In [3], it is shown that for complex Gaussian random variables zj and Z 2 
with identical variances and arbitrary means and covariance, the pairwise 
probability of error PrjlzJ 2 > |z,| 2 ]is given by 


p r{|z 2 | 2> |z,| 

2 } = |[l - QiJbM) + G(V^,V*)] 

where Q(a,P) is Marcum’s Q-function [4] and 


r " 

S x + S 2 - 2|p|V5^ cos(0, - 0 2 + 0) ^ S,-S 2 

W 2N.\ 

i-lpf V'-lpf. 


with 


(A-l) 


(A-2) 



P °2^( Z, " Z ‘)'( Zj ' Z2 ) 

0=arg{p}; 6 X = arg{zj}; 0 2 = arg{g (A . 3) 

Here we associate z^and Z 2 with w(A</>) and w(A<j>), respectively, of (23). 


Using (1) and (2) in (A-3), we get 

N-i-2 

_ N-l 

z, = V2PX^*-‘e ~° = 'j2Pe j * k ~ N + i N 


N-i-2 

__ N-l 

z 2 = V2 ?£«**-*« ~° = i j2Pe 5 

i«0 


(A-4) 


where 


RAGE 



INTENTION AfeLY BLANK 


PRECEDING PAGE BLANK NOT FILMED 
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N-i- 2 

A N - 1 j y y.( A * k ~‘- m ~ A *k-i-m) 

S=^e m=0 

1=0 

where it is understood that the summation equals zero if the upper 
summation index is negative. Substituting (A-4) into (A-3), gives 


(A-5) 


S l = PN 2 ; S 2 = P]<S| 2 
Also, using (3) 


(A-6) 


N= — E 
* 2 


r 

N-i-2 

r 


n-i ~i S ***-*-" 


< 

-° 



i=0 


V 


4 




(A-7) 


and 

f N-i-2 \ 

i yv-iiV-i >1 

^ iy z 1=0 n =0 

^ r N (A-8) 

Finally, substituting (A-6) - (A-8) into (A-l) gives the desired result, namely. 


Pr{f) > T]\A0} = Ul - Q(y[b, Va )- + Q(sfc, V* )] 


(A-9) 


where 


where E s = PT is the energy per data symbol. 


(A-10) 
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Appendix B 

Proof of Eqs. (42) and (48) 

Starting with the definition of 5 in (27), we now write it in the form 
5 = l + 5V ai (B-l) 

i= I 

where 

(B-2) 

11=1 

Thus, 5 is the sum of N unit vectors the first of whose arguments is zero and 
the rest of whose arguments are increasingly larger sums of the phase errors 
in accordance with (B-2). Note that the values of the accumulated phase 
errors, oq's, also range over the set ± 27im/M; m = 0 ,1, ..., M/2 - 1. We are 

interested in determining the various possible solutions for the S^'s such 
that the maximum value of the magnitude of 5, namely 1 5 1 max of (49) is 

achieved. 

For arbitrary M, there are four situations that achieve 1 5 1 max - We shall 
refer to these as cases 1, 2, 3, and 4 which are described as follows. 

Case 1 : All N-l vectors e) a i; i = 1, 2, ..., N-l must be collinear and equal to 
gj2jt/M Thus, oq = 2it/M; i = 1, 2, ..., N-l which, in accordance with (B-2), 
has the single solution 

= = fti = 3.4 ,.... N (B-3) 

Case 2 : All N-l vectors el a i; i = 1, 2, ..., N-l must be collinear and equal to 
e -j2it/M Thus, cq = - 2n/M; i = 1, 2, ..., N-l which, in accordance with (B-2), 
has the single solution 


fyk-N+2 — 


fyk-N+i = 0; j = 3, 4,..., A/ 


(B-4) 
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Case 3 : Any N-2 vectors d a i must be collinear and equal to e)° = 1 and the 

remaining vector must be equal to For this case there are N-l 

different solutions. For example, suppose first that aj = 2 tc/M, and cq = 0; 

i = 2, 3, N-l. Then, 


=~; =0;(= 4.5 N 


(B-5) 


Next, let a 2 = 2n/M, and cq = 0; i = 1, 3, 4, ..., N-l. Then, 


&Pk-N+ 3~ fyk-N+4— ^ ' fyk-N+i ~ 0> * ~ 2,5,6,. ..,N 

In general for = 27i/M, £=1,2, ..., N-2 and cq = 0; i = 1, 2, 
have the solution 


(B-6) 

..., N-l; i * 6, we 


o. _ 2n . su - 2n ■ 

°Yk-N+t + 1 — m » UYk-N+l+2 ~ ’ 


fyk-N+i = 0; I = 1,2,...£,£ + 3 


Finally, for oqq.j = 2n/M and cq = 0; i = 1, 2, ..., N-2, the solution is 


(B-7) 


ft i= 1,2 N-i 


(B-8) 


Case 4 : Any N-2 vectors e) a i must be collinear and equal to e^ = 1 and the 
remaining vector must be equal to e~) 2K /M p or this case there are again N-l 
different solutions which are identical to those described by (B-5) to (B-8) with 
2n/M replaced by -2 ji/M and vice versa. 

We note that for cases 3 and 4, N-2 of the solutions are characterized by 
having one &J> = 2k/M, one 8<|> = - 2n/M, and the rest of the 8<j>'s equal to zero. 
The remaining solution has one 8<J) = 2 ji/M and the rest of the 5<t>'s equal to 
zero. 


To compute the accumulated Hamming distance of (47), where w(u,u) 
corresponds only to those error sequences that result in 1 5 1 ma x/ we proceed 
as follows. We assume a Gray code bit to symbol assignment where 
&j> = ± 2 ji/M corresponds to an adjacent phase symbol error and thus a single 
bit error or a Hamming distance equal to 1. Also, a value 5q> = 0 implies no 
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symbol error or a Hamming distance equal to zero. Thus, the following 
accumulated Hamming distances occur for each of the four cases. 

Casel: 


w(u,u ) = 1 


(B-9) 


Case 2: 


w(u,u) = 1 


(B-10) 


Case 3: 


w(u,u) = (2 )(N - 2) + (1)(1) = 2{N - 2) + 1 


(B-ll) 


Case 4: 


w >(u,u) = (2 )(N - 2) + (1)(1) = 2(N - 2) + 1 (B-12) 

Finally, the accumulated Hamming distance is obtained by summing (B-9) 
through (B-12) which yields (for N > 2) 

£ w(u,u) = 4(N-l ) (B-13) 

which agrees with (48). 

For N = 2, Cases 3 and 4 do not occur since N - 2 = 0. Thus, the 
accumulated Hamming distance is merely the sum of (B-9) and (B-10) which 
yields 

£w(m,m) = 2 (B-14) 

A£-A£ 

in agreement with (47). 

For M=2, = el 71 = an( j thus cases 1 and 2 are one and the 

same and similarly for cases 3 and 4. Thus, for binary multiple bit DPSK, we 
have only half the solutions in which case (B-13) becomes 
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***** 


(B-15) 


which agrees with (42). Q.E.D. 
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Appendix C 

Evaluation of the Chemoff Bound on 
Pairwise Error Probability 

Here we derive an upper Chemoff bound on the pair-wise error 
probability Pr|^ r\ t > j. In particular, from the definition of the 

Chemoff bound [11] 



(C-l) 

where "E" denotes the expectation operator and Xq is the Chernoff parameter 
to be optimized. The remainder of this appendix is devoted to an evaluation 
of E:|e Xo( ^~ ,|t) [d0 (O |. For simplicity of notation, we shall drop the subscript "i" 

with the understanding that we are referring to the ith branch in the trellis 
paths under consideration. 

From the definitions of T} , and fj we can write the difference f) - 77 in the 
matrix form 


77 — 77 = vv* t Fvv 

where 

H A£)l n 0] 
w( A0) ' ~L° *1J 

and the "T" denotes the transpose operation. Here 

N-n - 2 

n-i ~i 

w{A$) = X r *-« e *"° 

*=0 

N-i i -2 

W-l 

w(d£) = £/•*_„<? "=° 

hmO 



(C-2) 

(C-3) 

(C-4) 


58 


with 


r k = j2Pe J{9 ** k) + n t 


(C-5) 


Substituting (C-2) into (C-l), we get an upper bound of the form given in 
[8: Eq. (21)] which can thus be evaluated as (see [12] for a derivation of the 
result) 

l — J • det(I-2A,£F) 
where I is the identity matrix. 


(C-6) 


^ = £'{hH A<p} = 


[e{w(A4)W}\ 


= V2 P\ 




N-n-2 

~J Afk-n-m 


m= 0 

n=0 

N-n-2 

N - 1 “ y 

m=0 

n=0 


and 


R = Ie{(vv-/£) * (w -/i) T l A^} = i 


2^ 2N 0 g * 


2^5 21VgiV 


T 

J * j 


(C-7) 


(C-8) 


with 


N-n-2 


N-n-2 


A N-l 1 ^( A *k-n-m~ a *k-n-m) ^ w _, / ^S+k-n-m 

8=^e mm0 =^e -° (C-9) 

*=*0 *=0 

To evaluate (C-6) we need to compute the determinant and inverse of 
(I- 2X 0 R*F). Using (C-3) and (C-8) these are evaluated, respectively, as 


n 2 n 2 

det(I-2A^F) = l-4A 2 0 ^ + 4A 2 0 ^ 


Nfc? 


(C-10) 


and 


(I-2A 0 R*F) _1 = 


det(I-2A^*F) 


. + nM -txM 

* f ^ f 


2a„m: 


1-2A, 




'• J 


(C-ll) 
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Finally, using (C-10) and (C-ll) in (C-6) gives 

E, A(l -^)[jV 2 -l^ 2 ] 

E^e Xo -* TF -\A^ = 
where E s = PT S and 



(C-12) 


A=2A„-^ (C-13) 

Reinserting the subscript "i" on 5 to denote the fact that (C-12) applies to the 
if h trellis branch, and substituting the result in (C-l) gives the desired 
Chernoff bound on pairwise error probability, namely. 



The value of X that minimizes the bound of (C-l 4) must be found by 
numerical evaluation. 
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Appendix D 

Asymptotic Evaluation of Certain 
Integrals of Analytic Functions 

Let f(z) and g(z) be meromorphic functions of z on the complex plane, C. 
Consider the integral 

/ r (y) = J g(z)exp(yy/(z))dz (D-l) 

r 

where the contour T does not pass through any singularity of f(z) or g(z) and 
we assume that y > 0. Although the contour is usually specified, we can 
clearly choose another contour Tq such that the integral in (D-l) remains 
unchanged, i.e., Ir(y) = Ipg(y) provided that the region bounded by T and Tq 

does not contain any of the singular points of the integrand. In particular, the 
selection of the contour Tq will be an important consideration in what 

follows. 

We now describe a procedure for evaluating (D-l) in the limit of large y, 
that is, we seek the first term of an asymptotic expansion of Ip(y) as y — » ©°. 
The following assumptions, which are applicable to the case at hand, are 
made: 

(i) The contour of integration Tq may be chosen to pass through at least one 
critical point zq of f(z), i.e., a point in C where f'fzg) = 0 and the prime denotes 
differentiation. The critical points of f(z) on Tq are assumed to be non- 
degenerate, i.e., f"(zo) * 0 if f'(zg) = 0. 

(ii) Im{f(zg)} is an absolute minimum of Im{f(z)) for z e To* 

(iii) The functions g(z), g(z)/f'(z), and g’(z)/f'(z) decay sufficiently rapidly as to 
ensure absolute convergence of the integrals that appear. 

Let {zj,* i e Til denote the set of critical points of f(z) that lie along Tq. Let 
X(z) be an infinitely differentiable function which is identically equal to unity 
in a small interval J j containing z\ and vanishes outside Jj z> ] i where Ji is a 
small neighborhood of z\ in Tq. Writing g(z) in the form 
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8( z ) = X( z )8( z ) + (l~X( z ))8( z ) (D-2) 

then we can easily see that Ip(y) can be written as a sum of integrals of the 
form 


A(y)= Jg(z)exp(./y/(z))rfz 


f 2 (7) = J g(z)exp(yy/(z))f/z 

h 


(D-3) 


where Ji contains a unique critical point of f(z), namely Zj, and Kj is a contour 
(possibly infinite) which contains no critical points of f(z). If Kj is unbounded, 
then the integral over Kj is understood as an improper integral in the 
standard manner. The contour Tq is composed of the union of 0j / Kj; i e q) . 
Also, the function g(z) vanishes at the end points of Ji and Kj. We now 
investigate the asymptotic behavior of Ij(y) and l2(y). 


To evaluate 1 1 (y), we should map the contour line integral into an integral 
along a portion of the real line parameterized by the variable t. Thus, let £-(t) 

be a continuously differentiable mapping, with nowhere vanishing 
derivative, of the interval t e (0,1) onto the contour Kj. Then, keeping in 
mind that Kj contains no critical points of f(z), the first integral in (D-3) can be 
written as 


AM = (0)e*p(jr/(f «»)f ' (Ddt 


.^MpO'r /<C(»») 


*({«))* 


l f' i(CM) 
;> J ° /'(£(<)> 


dexp(;y/(£(0)) 


Integrating by parts results in 


(D-4) 


A(y) 


g(C ( f )) ex p(yy/(C (0)) 


Jr /'((«)) 




( g(C(Q) 


dt 


(D-5) 
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The first term in (D-5) vanishes in view of the assumption of rapid decay of 
g(z)/f'(z) if Kj is infinite, and the vanishing of g(z) at the end points t = 0,1. 

Therefore, 

(D-6) 

whose absolute value is upper bounded by 

IwrH * ^J>p(-rin.{/(f (D-7) 

Finally, since by assumption (ii), 

minlm{/(z)} = Im{/(z i )}=c i (D-8) 

then, in the integration interval of (D-7), 

exp(-yIm{/(C(0)}) £ exp(-yc,) (D-9) 

and hence. 


|A( y )| ^ C, y exp(- yc,); C, = 



(D-10) 


The second integral of (D-3) involves evaluation in the neighborhood of 
the critical point z[ under the assumption that the critical point z\ is a 
minimum of Im{f(z)}, z e Tq. The asymptotic expansion of this integral is 
given in [13: Chap. 8]. Using the first term as an approximation to 12(7), we 
obtain 

7 2 (y) = ^ yf"{ z ) eX P 0'y/( Z <>)g( Z «> 

Once again considering the absolute value of this integral, we get 


Uy] = d^cM-yc,); D i = 


2 nj 




■«(M 


(D-12) 


which clearly dominates I Ii(y) I of (D-10). It is easy to see that the 
contribution of the critical points which are not minimum of Im{f(z)} are 
dominated by 1 12(7) I • Therefore, for large y, Ip(y) is asymptotically given by 
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A- ( 7 ) = Z Jt^TT exp(yy/(z i ))g(z i ) (D- 13 ) 

where the summation is over all the critical points z[ of f(z) that are absolute 
minima of Im{f(z)} on Tg. Note that if z g is the only such critical point along 
the contour Tg, then the summation in (D-13) contains only a single term. 

As we shall see, in special cases of the above general theory, it is possible to 
choose the contour To such that this is true and hence the numerical 
evaluation is considerably simplified. 

One very special case of the above occurs when f(z) and g(z) are of the form 
[see (73)] 


/(*)=£- 


k + l (D-14) 

where ^i. Pi, and P2i are real coefficients (note that = 1) and the contour T 
is the straight line extending from z = - °° + je to z = + je. It is sufficient to 

have e > 0 to avoid the pole of g(z) at the origin. However, as mentioned 
above we should redefine the contour (to Tg) so that it passes through a 
critical point of f(z). For B = 1, it can be shown that there are two critical 
points both of which lie along the imaginary z axis, i.e., zj = jpj; i = 1,2. One 
of these critical points, say pi, lies below the singularity of g(z) at z = = j 

whereas the second critical point, say P 2 , lies above this singularity. Thus, in 
accordance with the above, if we choose e = pi then the contour Tg will pass 
through only one critical point of f(z) and not cross any of the singularities. 
Hence, the summation in (D-13) reduces to just a single term. The details of 
the evaluation are carried out in Section 8.2 of the main body of the report. 

When B > 1, then 

/' ( Z ) = -j\l ~ 11M? Z hi = g 

1 ki (S¥+if 

(D-15) 

which, in general, yields 2 + 4(D-1) critical points of f(z). Since it is not, in 
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general, true that all the critical points lie on the imaginary z axis, one must 
evaluate these critical points numerically for each particular case and then 
select the contour Tg appropriately as discussed above. It can be shown 
however (see Appendix E), that for arbitrary B, there exists a unique critical 
point, say zg, on the imaginary z axis in the interval 0 < z < j. Therefore, 
since <, 1 for 1 <; i < B (see the definition of in Eq. (72) of the main text), 
and thus, aside from the singularity of g(z) at z = 0, all of the singularities of 
g(z) and f(z) occur along the imaginary z axis in the interval j < z < j°o, it is 
sufficient to find this unique critical point zg = jyg (0 < yg < 1) and choose the 
contour Tg to pass through it. . . . 

It remains to satisfy condition (ii). Let v(z) = Im{f(z)}. Then, from the 
computations in Appendix E, it is clear that 


d\(x+Jy) 

dy 2 


<0 


(D-16) 


for 0 < y < 1. Since v(z) is harmonic, i.e., for all analytic functions f(z), we 
have 


d*v(x + jy) , dMx + jy) _ n 

dy 2 dx 2 

then 


(D-17) 


dMx±Jy) 

dx 2 


>0 


(D-18) 


and thus the critical point zg = jyg is a local minimum of v(z) along the line 
z = x + jyg parallel to the x axis. It is straightforward to show that the contour 
Tg can be deformed to a contour Tg' such that 


IrSY) = I r (Y) (D-19) 

* *0 

and where zg = jyg is the global minumum of v(z) on Tg'. Thus, condition 
(ii) is satisfied for the critical point zg. 
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Appendix E 

Proof That There Exists a Unique Critical Point 
of f(z) Along the Imaginary z Axis in the 
Interval 0 < z < j 

Here we present a proof that the function f(z) (defined in Eq. (73) of the 
main text or Eq. (D-14) of Appendix D) has a unique critical point (i.e., the 
value of z where f’(z) = 0) along the imaginary z axis in the interval 0 < z < j 
Consider writing the derivative f'(z) along the imaginary z axis in the form 

= (E-i) 

i=i 

where, from (D-15), 


f'.(jy) = E.Pl&l W& y + @ 2! 

‘ (i -tyf 

Evaluating (E-2) at y = 0, we have 




At y = 1, (E-2) becomes 


n 2 -\51 2 




(E-2) 


(E-3) 



(E-4) 

The first factor in (E-4) is clearly positive. We now wish to show that the 
second (bracketed) factor in (E-4) is negative. The steps are as follows: 
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(E-5) 

Thus, since each term in f'(jy) changes sign in the interval 0 < y < 1, there 
must be at least one value of y at which f'(jy) = 0, i.e., at least one critical point 
of f(z) occurs in the interval 0 < z < j. It now remains to show that there is 
only one such point, i.e., the critical point of f(z) that occurs in the interval 
0 < z < j is unique. We shall do this by showing that f"(z) is always negative 
in the interval 0 <.z<) and thus f’(z) is monotonically decreasing in this same 
interval. 

Consider the derivative of (E-2). After considerable simplification and use 
of the definitions of Pj, P2h and £i in (72), we obtain 



(E-6) 

It is straightforward to show that the sum of the first, third, and fifth terms of 
the numerator are non-negative in the interval 0 ^ v < 1 (or equivalently 
0 < y < (£i)”l)- In particular, 
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-i(n 7 - |SJ 2 ) ,/2 v 5 - 4(tt 2 - |5,f) 1/2 v 3 + 6(jV 2 - |<5 ,| 2 ) 1/2 v 
= -2(iV 2 - |5j 2 ) 1/2 v(v 4 + 2v 2 - 3) > 0 for 0 < v < 1 

Thus, 


-2 (N 1 - |<5,f f V - 4(W ! - j<5f) 1/2 v 3 + «(n 2 ■ - |<S,f )' ,! v 

£ —2Nv 5 - 4Nv 3 + 6Nv forO<v<l 
Substituting (E-8) in (E-6), we get 


(E-8) 


fm* 


-2Nv 5 + 6Nv 4 - 4Nv 3 - 4 Nv 2 + 6 Nv - 2 N 
(l-vt 


_ -2 A(v 5 - 3v 4 + 2v 3 + 2v 2 - 3v + 1) 

(i-v 1 )* 


—2N(l - v) 3 

■hT 


SO; 


0£v£l 


(E-9) 


Finally, since the derivative of (E-l) is the sum of the derivatives of each term 
in the summation, and recalling from (72) that ^ < 1, then from (E-9) we 
obtain the desired result, namely. 


/;(/y)<0; 0<y <1 (E-10) 

or, equivalently, fi'(jy) is monotonically decreasing in 0 < y < 1. Q.E.D. 


Although not of specific interest, it can also be shown that the critical 
points of f(z) that do not lie on the imaginary z axis occur in pairs that are 
symmetric about this axis. That is, if zq = xg + jyo is a critical point, i.e., f'(zg) = 
0, then z\ = -xQ+jyo satisfies f’fzj) = 0. 
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